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We use nonlinear perturbation theory to obtain new solutions for brane world models that incor- 
porate a massive bulk scalar field. We then consider tensor perturbations and show that Newtonian 
gravity is recovered on the brane for both a light scalar field and for a bulk field with large negative 
mass. This latter result points to the viability of higher-derivative theories of gravity in the context 
of bulk extra dimensions. 



I. INTRODUCTION 

One of the unsatisfying features of the Randall- 
Sundrum models has always been the need to fine-tune 
the brane tension and the bulk cosmological constant. 
Immediately after these models were proposed, many 
authors considered adding bulk scalar fields in order 
to evade this problem. It was hoped that the brane- 
bulk system would also be stabilized by the presence of 
these fields. However, it has proved difficult to find self- 
consistent analytic solutions to all but the simplest sce- 
narios, e.g. superpotential models P, . Furthermore, 
the order of complexity rises when one considers a fully 
time-dependent scenario. For this reason, work has re- 
cently been concentrated on numerical simulations 4]. 

On the other hand, the insight of Randall and Sundrum 
was that gravity can still be effectively 4-dimensional de- 
spite the presence of a fifth, spatial dimension. There- 
fore, if one complicates their simple models by including 
a bulk scalar field, it is vital to ensure that gravity has 
the usual Newtonian limit, at the very least. 

In this paper, we tackle analytically the case of brane 
worlds with a massive bulk scalar field. We conduct two 
distinct perturbative analyses. In the first, we utilize 
ideas in nonlinear perturbation theory to obtain leading 
corrections to the background warp factor and bulk field; 
it is the nature of the former that is crucial to the local- 
ization of gravity on the brane. Then, we consider metric 
and matter perturbations of the brane-bulk system, and 
compute the Green's function in response to a point-mass 
on the brane. We show that the Green's function has the 
required 1/r-behavior in the long-distance limit. 

Although our presentation here is something of a math- 
ematical proof-of-concept, it sets the stage for a subse- 
quent paper that will explore these issues in the context 
of theories of higher-derivative gravity. 



II. BACKGROUND SOLUTIONS 

We begin our analysis by considering the model of one 
Minkowski brane in a bulk that contains a massive scalar 
field <j). The scalar field potential is V = \M 2 cjP. The 
brane is coupled to the scalar field via its tension; to be 
specific, we choose a Liouville potential for the tension, 
X{4>) = e - 2 */^3. I n fact, the explicit choice of \(<f>) is 
not so important since we only ever need its value at 



0(0). Our choice effectively fixes the degree of freedom 
associated with X'(4>(0)). 

The line-element for the brane-bulk system may be 
written as 

ds 2 = e 2X ^r) liV dx' i dx v - dy 2 , (1) 

with the brane at y = 0. The equations of motion that 
follow from the Einstein field equations are 

-3X yy -6X 2 = \4> 2 V + \M 2 ^ 2 , (2) 

6X 2 = \4>l ~ \M 2 4>\ (3) 

(pyy + AX y (f> y = M 2 (j). (4) 

Note that we will often use a subscript to denote a deriva- 
tive when there is no cause for confusion. The junction 
conditions at the brane position are 

X y (0) = -±, ^(0) = -^=, (5) 

where we have assumed Z2-symmetry — we evaluate 
quantities as y J, — and used the fact that 0(0) = 0, 
which follows from eq. (pj- Without loss of generality, 
we choose A"(0) = 0. Finally, we point out that we have 
made all quantities dimensionless by rescaling in terms 
of a characteristic length A _1 (0)Mp and characteristic 

scalar field value Mp 2 , where Mp is the 5-dimensional 
Planck mass. 

A. Small M 2 

We let e = M 2 be our perturbative parameter. Then, 
we introduce a "strained" coordinate, 

x = (l + f ie)- 1 y, (6) 

in place of y. This technique, known as the Lindstedt- 
Poincare method can make a divergent perturbative 
expansion convergent by an appropriate choice of fi. 

We now let X = X (x) + eXi(x) + ... and — 4>q{x) + 
t<f>i(x) + . . ., and expand the equations of motion and 
junction conditions up to first order in e. The zeroth order 
equations yield 

X = i In (1 - fx) , O = 2V3X . (7) 

The obvious feature of the solution is the singularity at 
x = | . It is easy to check that this is not a coordinate sin- 
gularity; indeed, such singularities appear to be generic 
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in scalar field models Q, Q • We will discuss what to do 
about them in subsection III CI 

The first order equations are much more involved and 
we omit the expressions for X\ and <t>i- We are not pri- 
marily interested in their contribution to the overall so- 
lution at this stage, but rather in what they imply about 
/i. To make this last point more transparent, it is enough 
to know that the homogeneous solution to X x is 

X[ h) = AO*) + M- t (8) 

where A and B are constants. The problem is that the 
latter term will dominate X : even for small e, the warp 
factor e x will pick up a term like (1 — |x)~ 3 / 4 . Such a 
term is sometimes called a secular term. The solution is 
to choose so that B{^i) = 0. It turns out that 

M = -i (9) 

and then 

Note that, even at zeroth order, an effect of the perturba- 
tion is seen. However, we now need the first order terms 
to ensure the junction conditions are satisfied. That the 
singularity occurs at y = |(1 — ^e), is borne out by nu- 
merical simulations. 



lower order in e. The solutions of the equations resulting 
from substituting the above expansions into the equa- 
tions of motion and the junction conditions are 

Ei(u) — sin 2 u, Fi(u) = -4= sinu. (13) 



C. Eliminating the singularity 

As pointed out above, the presence of singularities in 
our solutions is both generic and problematic. The stan- 
dard method of dealing with them is to add a "regulator" 
brane to the system. If the singularity is at y = y s , then 
placing the second brane at y — y+ < y s will, due to Z2- 
symmetry, "slice off" the section y > y+, which includes 
the singularity. 

The insertion point y* is not arbitrary, however. If we 
let the brane tension of the regulator brane be A* (</>), 
then a second set of junction conditions which must be 
satisfied are 

xM = \x i< {ct>), ct>M = -\Ki<t>), (14) 

where all quantities are evaluated as y | y*. In what 
follows, we will assume that these conditions are met. 
The subsequent question of whether the two-brane setup 
can be stabilized is left for a future numerical study. 

III. TENSOR PERTURBATIONS 



B. Large —M 2 

In this case, the convenient perturbative parameter is 
e = The Lindstedt-Poincare method is not re- 

quired here, but because the e — > limit is problematic, 
it is also not obvious what perturbative expansion will 
work. From numerical observations, it is apparent that 
X ~ \ ln(l — |y) but that highly oscillatory terms arise 
in X y and X yy . Curiously, the solution for X is what one 
would obtain in the usual Randall-Sundrum case with 
no bulk cosmological constant if eq. © could be ig- 
nored. Therefore, the limit we are considering here can 
be thought of as the best a scalar field can do to mimic 
A = 0. Note that, once again, there is a singularity at 
finite coordinate distance from the brane. 

Expansions for X and <j) that capture the oscillatory 
nature of the solutions are as follows: 

00 

X = lln(l-ly)+J2e^(l-ly)- 2n E n ^),(ll) 
n=l 

00 

- E^a-^r 2 "^®- (12) 

n=l 

In fact, for our purposes, the first term under the sum- 
mation will be sufficient. The important point, however, 
is that taking derivatives of these terms makes them of 



Having obtained the background solutions, we are now 
interested in perturbing about them. In particular, we 
focus on tensor perturbations in the brane-bulk system 
in order to understand the nature of gravity for a brane- 
bound observer. 

The most convenient framework in which to describe 
perturbations is the explicitly conformal line-element 

ds 2 = a 2 (Y) (n^dx^dx" - dY 2 ) (15) 

that follows from transforming line-element (I) to a 
new coordinate: Y = f e~ x ^ v ^dy'. The warp factor is 
a(Y) — e x ( y \ The gauge-invariant tensor perturbations 
decouple from the scalar and vector perturbations, and 
take the form 

d 9 P=a 2 h^, 6g$=8gjjP = 0, (16) 

where indices on are raised (lowered) by tj^ (r]^) 
and h^u is transverse traceless: 

h% = dphTu = 0. (17) 

The only wrinkle in this otherwise straightforward 
setup is that one must also consider perturbation in the 
brane positions @|. To be precise, there is a scalar de- 
gree of freedom associated with each brane, which adds 
an additional term to the perturbed junction conditions. 
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If there were only one brane and no matter perturba- 
tion on it, the perturbation in position could simply be 
eliminated by a redefinition of the coordinates. In our 
scenario, we will suppose there is no matter perturbation 
on the second brane and consider its position fixed. On 
the other hand, we wish to calculate the gravitational 
response to a point-mass on the first brane and thus its 
position, here corresponding to the radion, must be per- 
turbed. 



A. Equations of motion 

In terms of our new coordinate, the original unper- 
turbed brane is located Y = 0. We now suppose it is 
perturbed to Y = Ci x>1 ) — to first order, y — C,(x^) also — 
following a general matter perturbation given by the 
energy- momentum tensor t^ v . The resulting equations 
of motion for the tensor modes are 

(d 2 + 3^-d Y - □) V = 0, (18) 

subject to the junction conditions 

dyhfiulyio = T tu> ~ h^iivT - = o>„, (19) 

Because there are two conditions, we expect only a dis- 
crete set of gravitational modes. 

The fact that the left-hand side of l|19JI is traceless 
gives us the important equation for the displacement of 
the brane @, namely 

□C = -±r. (21) 

If we assume d^r^v = 0, then the transverse nature of 
/i M „ is reflected in eq. (12H as well. In principle, we must 
solve this equation and substitute the solution into eq. 
(tlflll . However, in this paper, we will content ourselves 
with computing the appropriate Green's function from 
which can be constructed. 

B. Green's function solution 

Equations iflSjl , iflfljl and ll2f))) can be conveniently com- 
bined as 

(d Y + 3^-d Y - □) V = 2a^(x)6(Y), (22) 

allowing us to write a formal solution in terms of the 
retarded Green's function given by 

(d Y + 3—d Y - □) G(x, x 1 , Y, Y') = 5 {i) {x-x')8{Y-Y'). 

(23) 

After a Fourier transform, this becomes a 1-dimensional 
Green's function problem 

D k G(Y, Y') = (d Y + 3^-d Y + k 2 ^j G = 5{Y-Y'), (24) 



which can be solved in terms of the complete set of nor- 
malized eigenfunctions of D k . 

To be precise, we suppose D k Xm = (k 2 — m 2 )xm, with 
boundary conditions <9yx™(0) = d Y Xm(Y-,<) = 0. (These 
boundary conditions reflect both the continuity of d Y \m 
and the requirement of ^-symmetry.) Then, introducing 
"0m = s 3 ' 2 Xm, we obtain the Schrodinger-like equation 

(d Y+m 2 -^^^ m = 0. (25) 

This is usefully rewritten as Q 

m 2 ) V™ = 0, (26) 

where 

D ± = 8 Y ± \X Y . (27) 

The boundary conditions on the tp m can be written as 
Z)_V> m (0) = D^ m (Y i ,) = 0, and the eigenfunction nor- 
malization is 2 fy* dYip m tp n — S mn . Then, 

G(YY>) = Ya- ^ Y 2 ) ^P . (28) 

m 

The information we particularly want from the Green's 
function is the response, in the stationary limit, at posi- 
tion x on the brane due to a disturbance at x' also on 
the brane. We have 



G(x,x') = 




= -i^E^(0)e—, (29) 
m 

where r = |x — x'|. A similar analysis was carried out in 

C. Recovering the Newtonian limit 

Inspection of eq. (129)) shows we obtain the necessary 
leading 1/r-behavior if there is a zero mode. The contri- 
bution of other modes may spoil this result if there are 
tachyon modes, or if there are modes with mr <g; 1. This 
latter scenario is model-dependent but we can show Q 
that there is always a zero mode and never a tachyon 
mode, i.e. m 2 > 0. 

Multiplying eq. l(2fij) by ip m and integrating over Y 
gives 

m 2 = -2 f dYip m D + D-il> m 
Jo 

= -2i\} m D-ii m +2 / dY(D_iP m ) 2 . (30) 
Jo 
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The first term is zero by the boundary conditions and 
the second term is non-negative. Furthermore, the zero 
mode is given by D-tyo — 0, that is 



V>o oc a 3/2 . 



(31) 



IV. 



GRAVITY AND A MASSIVE BULK 
SCALAR FIELD 



We now apply the above formalism to the specific back- 
ground solutions computed in section [H] In both cases, 
small M 2 and large — M 2 , we will concentrate on the sit- 
uation at zeroth order in e. For a start, as we have seen, 
nonlinear perturbation theory already gives perturbative 
effects at zeroth order, and secondly, because eq. 12511 is 
of Schrodinger-type, it is straightforward, though hardly 
illuminating, to compute perturbative corrections to the 
if) m . We will assume that Y+ -C r and call this the long- 
distance limit. 



A. Small M 

In this limit, we cannot completely ignore the first or- 
der term in our solution for X because it is necessary for 
getting the boundary condition at Y = right. How- 
ever, to a good approximation, we may account for X\ 
by simply rescaling Xq : 



X ~ i£Mn (l - fB^y) , 



where B = 1 — g£. It follows that 



Y s /6 



(32) 



(33) 



with Y s 



\e to first order in e. The solutions to eq. 



H25I1 are then built from linear combinations of Bessel 
functions of the first and second kind: 



V> ro = (Y a - Yf' 2 J p (m(Y s - Y)) , 



(34) 



where p = j(Y s — 2) ~ ~| e - The boundary conditions 
become J p +i{mY s ) = ^ p +i(m(Y s — Y*)) = 0. In the limit 
of large mY s , we find 



Am = irY* 1 > ttY s 1 > r 1 , 



(35) 



and therefore that the infinite tower of states will not 
seriously alter ll29")l . The Newtonian limit will also not 
be affected by light modes: there is no massive bound 
state in the limit of small mY s , and we conclude that the 
lightest mode always satisfies 



mi ~ y: 



(36) 



B. Large —M 2 

We can use the zeroth order term for X to com- 
pute the coordinate Y. However, the second order term 
E\(y(Y)/e) is crucial to the computation of the effec- 
tive potential in i!25|) : after being differentiated twice, it 
becomes a zeroth order quantity. Actually, it is slightly 
more convenient to use the physical coordinate y. Letting 
Cm = a 1/2 ip m , we have 



cos 



2.u 



1 



(i-!) 2 



Cm = 0. 



(37) 



There is some hope of solving this equation in the limit 
that the regulator brane is brought very close to the first 
brane. To be specific, we suppose ~ e. Then letting 
y = ez, we obtain a Mathieu equation 



(d 2 + A- 2q cos 2z) Cm = 0, 



(38) 



where A — e 2 m 2 and q = —ike 2 , subject to d z C 



2 | > ..•■>■ 'J- ■■ ■ sm (0) 

If we could find bound states with m ~ 



^zCm(^-k) 3 

0(e n ), n > — 1, it might be that raj/ t <C mr <C 1, and this 
would be fatal for the Newtonian limit. As it happens, 
this does not occur. The solutions @ to © are linear 
combinations of F v (±z), which for small q become 



2iz 



-Hz 



(39) 



where A ~ v 2 + \q 2 / (v 2 — 1) for v real and non-integer. 
Then, one can show that it is not possible to satisfy the 
boundary conditions for small, real v. This implies that 
either v ~ 1 or v is purely imaginary, i.e., in an instability 
band of the Mathieu equation. In either case, A ~ 1. 



V. CONCLUSIONS 

We have conducted two distinct perturbative analyses 
of brane world models that incorporate a massive bulk 
scalar field. First, we used nonlinear perturbation the- 
ory to find background solutions in the limits of small 
mass and large negative mass. In both cases, we found 
it was necessary to add a second brane to eliminate a 
singularity. Second, we considered tensor perturbations 
about these background solutions. We were able to show 
that the brane-bound observer will observe only slight 
modifications to Newtonian gravity in the case of a light 
bulk scalar field. This result also holds true for a bulk 
scalar field with large negative mass, which in some ways 
mimics a bulk cosmological constant. Even though such 
a bulk potential may be unphysical, it arises naturally in 
higher-derivative theories of gravity that are conformally 
related to scalar-tensor theories, including general rela- 
tivity with scalar field matter. Thus our analysis sheds 
light on the viability of these theories in the context of 
bulk extra dimensions. We will make this connection 
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more rigorous in upcoming work We also look to SFB 375. 

study the stability of the two-brane system introduced 

here. 
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